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We discuss a correlated two-photon imaging apparatus that is capable of producing images that 
are free of the effects of odd-order aberration introduced by the optical system. We show that both 
quantum-entangled and classically correlated light sources are capable of producing the desired 
spatial-aberration cancellation. 
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I. INTRODUCTION 

Two-photon imaging, also known as ghost imaging, in- 
volves the use of coincidence measurements to form im- 
ages via photons that never interacted with the object 
being viewed. It has been a topic of great interest since 
its discovery using entangled photon pairs Initially, 
it was believed that the entanglement was a necessary 
ingredient for the effect, but it has since been found that 
most aspects of ghost imaging can be simulated using 
spatially-correlated classical light 0, 0] , including ther- 
mal and speckle sources 0-H|. 

Separately, it has been demonstrated that the entan- 
gled photon pairs produced in spontaneous parametric 
downconversion (SPDC) may also be used to cancel some 
of the effects of frequency dispersion flfjl - H2j | or spatial 
dispersion (aberration) [13l - [l5j in interferometer experi- 
ments. 

The goal here is to show that two-photon imaging can 
be done with spatially correlated pairs of light beams in 
such a way that odd-order aberration effects introduced 
by an imaging system may be cancelled. This may be 
done with either classically correlated beams or quan- 
tum entangled photon pairs. The key idea is to partially 
collapse the two arms of the ghost-imaging setup so that 
both of the output beams pass through the same optical 
system in an anticorrelated manner, but with only one 
of the beams passing through the object. Here we will 
illustrate the method using the simplest possible imaging 
system, a single lens. 

The outline of the paper is as follows. In section II, we 
review the basic idea of correlated two-photon imaging. 
In section III we discuss aberration-cancelled two-photon 
imaging with an entangled light source. In section IV we 
briefly discuss how to do the same with a classical source, 
followed by conclusions in section V. 



II. TWO-PHOTON GHOST IMAGING 

Two-photon correlated imaging, or ghost imaging 
is done with an apparatus like the one depicted schemat- 
ically in figure 1. In the original version, the correlated 
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FIG. I: (Color online) Schematic depiction of two-photon 
(ghost) imaging setup used to view object G. 



photon source is a x*- 2 - 1 nonlinear crystal pumped by a 
laser, leading to spontaneous parametric downconver- 
sion. This is the source we will use in section III. Entan- 
gled photon pairs with anticorrelated transverse spatial 
momenta q and — q travel along the two arms of the appa- 
ratus. The object to be viewed introduces a modulation 
(either the transmittance or reflectance of the object), 
given by a function G(x). This object is placed in arm 2 
(the upper branch), followed by a large bucket detector, 
D 2 . The detector's area is integrated over, so D 2 can not 
record any information on the position or momentum of 
the photon that reached the object; all this detector is 
able to tell us is whether the photon reached the detec- 
tor unimpeded, or whether its passage was blocked by the 
object. In the other branch of the apparatus there is no 
object, and all of the photons reach an array of pointlikc 
detectors or some other form of spatially resolving de- 
tector without hindrance. A coincidence circuit is used 
to record a count every time a photon detection occurs 
simultaneously (within the coincidence time window) at 
both detectors. By plotting the coincidence rate as a 
function of position x\ in detector f , we build up an im- 
age of the object. This is true even though photons that 
actually encountered the object in branch 2 left no record 
of the object's position, and the photons in branch f that 
do carry position information never encounter the object. 
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FIG. 2: (Color online) Schematic depiction of setup for odd- 
order aberration- cancelled imaging of object G. 



FIG. 3: (Color online) Implementation of scheme of fig. 2 
using parametric downconversion. 



The essential ingredient is the spatial correlation of 
the downconverted photon pair. The question arose as 
to whether the entanglement of the photons was neces- 
sary, or if a classical source with anticorrelated transverse 
momenta could mimic the effect. It was found 0, 9 that 
this was indeed possible. The correlated light source in 
this case consists of a beam steering modulator (a rotat- 
ing mirror, for example) which directs a classical light 
beam through a range of q vectors, illuminating different 
spots on the object. The beamsplitter then turns the sin- 
gle beam of transverse momentum q into a pair of beams 
with momenta q and — q. The results were similar to 
those with the entangled source, but with half the visi- 
bility. It has since been shown that spatial correlations 
present in radiation produced using thermal and speckle 
sources may also lead to ghost imaging (0-[j|). 



III. ODD-ORDER ABERRATION 
CANCELLATION IN TWO-PHOTON IMAGING 

Many optical devices work by adding position- 
dependent phase shifts to a beam as it passes through the 
device. A lens, for example, produces focusing by adding 

ik 2 

to the beam a quadratic phase shift e 2 f , where k, f 
and x are the wavenumber of the beam, the focal length 
of the lens, and the displacement of the given point in 
the beam from the axis of the lens. However, imperfec- 
tions in the shape of the lens and variations in the index 
of refraction of the material from which it is made may 
lead to additional unwanted phase shifts e Ic ^ x ' beyond 
those intended. These unwanted phases are an example 
of optical aberration: they lead to distortions of the out- 
going wavefronts, and consequently produce distortion 
in the final image. The aberration function </>(x) may be 
parameterized in a number of ways, for example by ex- 
panding in Zernike or Seidel polynomials [16Hl8j . Here 
we are uninterested in the details of how the function is 
represented and are concerned only with the fact that it 
may be split into a sum of parts which are either even or 



odd under reflection about the axis: 



where 



(x) = 4> even (x) + f> 0£i d(x), 



fieven ( x) — 4*even (x) j 
4>odd(~*-) = -</>odd(x). 



(1) 



(2) 
(3) 



The even-order terms include astigmatism and spheri- 
cal aberration, while coma contributes to the odd-order 
terms. 

We now wish to take the two-photon imaging setup of 
fig. 1 and alter it in order to cancel as much of the aber- 
ration of the lens as possible. Note that we should aim 
only to cancel the odd part of the aberration, since can- 
cellation of the even-order phases would also completely 
cancel the effect (the quadratic phase shift) of the lens 
itself. 

So consider the setup shown in fig. 2. This is similar 
to fig. 1, but with one main difference: we have partially 
merged the two arms, separating the photons in each pair 
only after they have passed through the lens. The two 
photons both arise from the same well-localized point £ in 
the source, but emerge with opposite transverse momenta 
±q. So both pass through the same lens, but tend to pass 
through it on opposite sides of the axis. 

Fig. [3] shows a slightly more detailed picture of a par- 
ticular embodiment of this scheme, using entangled pho- 
ton pairs produced via parametric downconversion. Sup- 
pose the field in the pump beam as it enters the crystal is 
E p (£). We arrange for the distances to satisfy the imag- 
ing condition 
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(4) 



The lens will have a pupil function which we assume to 
consist of the usual quadratic phase plus an additional 
phase function describing aberrations: 



(x) 



(5) 
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FIG. 4: (Color online) Unfolded version of the apparatus in 
fig. 3, with the signal and idler drawn as a single continuous 
ray. 

The beam splitter sends half of the photons to a spatially- 
resolving detector D\, the other half to a large bucket 
detector D 2 , with no spatial resolution. The object in- 
troduces a modulation G(x2) only in branch 2. The im- 
pulse response function for beam j (J = 1,2) to travel 
from transverse position £ in the crystal plane to trans- 
verse position Xj in the detection plane is then 

h(£,Xj) = J ff 3 (x j )e lfc [ (?2/Zl)+(X3/Z2) ] /2 (6) 

X f e -' fc [(«/ z i) + ( x j/^)]-x' e #(x') d 2 a ,/ 



where 



G(x 2 ) 



for j = 1 
for j = 2. 



(7) 



One further view of the apparatus is useful. Fig. @] 
shows an unfolded version of the setup in the Klyshko 
picture [l9|, [2(| . In this picture, the signal and idler can 
be viewed as a single continuous ray from one detector to 
the other. In the approximation of a plane wave pump, 
this ray is undeflected at the crystal because of the per- 
fect anti-correlation between the signal and the idler wave 
vector directions (q = — q), causing the signal and idler 
to hit the lens at a pair of points arranged symmetri- 
cally about the point where the ray crosses the crystal, 
x' - £ = £ - x". 

Assuming a thin crystal and narrow-band filters in the 
beams, the coincidence amplitude (or equivalently, the 
two-photon wavefunction) in the detection plane can be 
written (up to overall normalization) in the simple form 

in 

<A(xi,x 2 ) = j J B p (£)/ ll (£,x 1 )/i 2 (£,x 2 )d 2 £ (8) 

E p {£)G^2)e {lke/zi) e {tk,2z2)( y x * +x i) 

x e -ifc£-(x'+x")/zi e -ife(xi-x'+X2-x") 

x e l (^+^"^d 2 (d 2 x'd 2 x". (9) 



The coincidence rate is 

R{xx) = J d 2 x 2 |^(xi,x 2 )| 2 

If we define 

nz) = Ep{tt)e* e,z \ 

then the amplitude may be written in the form 
V(xi,x 2 ) = G(x 2 )e 2fe (^ 2 + X 1 2 )/ 2 ^ 
F 



(10) 



(11) 



(12) 



fc(x +X ) \ ife( xl . x ' +X2 . x ")/ Z2 



x e i (^ + ^)d 2 x'd 2 x"d 2 ^. 

For simplicity, we take the pump to be an approximate 
plane wave over the extent of the object, so that we may 
set E p (€) equal to a constant. We also assume that we 
are working in the far field, where the £ 2 term in the ex- 
ponential may be neglected in comparison with the other 

terms. In that case, we find F (x' + x")^ is propor- 
tional to 5^ (x' + x"), so that the amplitude reduces to 

iKxi.xa) - G(x 2 )e ifc (** 2 +*i 2 )/(2* 2 ) (13) 

e -ifc(xi-x 2 )-x' / z 2 e 2i0„ c „(x') ^2 x > 



(14) 



X $ 

where we have defined 



G(x 2 )e 4fc ( X22+Xl2 V( 2 ^ 

k (xi - x 2 ) 



-2 



= 2# c „ e „(x) 



(15) 



and the tilde denotes the Fourier transform. Given this, 
the coincidence rate is then 



R(xi) = J d 2 x 2 |V'(xi,x 2 ) 



d x 2 



G(x 2 )$ e „ 



• 11 1 (xi - x 2 ) 

z 2 



(16) 



(17) 



Eq. [T7] is our main result. Note that it depends only 
on the even-order aberrations; all odd-order terms have 
dropped out. In the special case that there is no aberra- 
tion in the lens, this reduces to 



i?( Xl ) = |G( Xl )r 



(18) 



It is apparent that we will always obtain unit magnifica- 
tion. 

The odd-order aberration cancellation is exact only in 
the far field and in the case of a plane wave pump. As 
the distances involved decrease or as the pump amplitude 

deviates from a constant, the factor F (j^ (x' + x")^ will 

no longer be a delta function, so that the aberration can- 
cellation will become only approximate. 



4 



Illumination 
^ 




Object 
G($) 



FIG. 5: (Color online) Simple incoherent imaging system. 
The distances z\ and Z2 obey the imaging condition — + — = 
l 
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The coincidence rate may be compared directly to 
the image intensity I(x\) collected by an incoherently- 
illuminated single-lens imaging system such as that of 
fig. 5. Up to overall normalization, the output inten- 
sity of the system in fig. 5, taking lens aberrations into 
account, is 



j( Xl ) = / d 2 t;\G(OE P (Of 



(19) 



x / d 2 x' ' e ~ ikx '-(KMHO*!/ 'e'^C*') 



~1 



xi 
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where 



$(x) 



J<f>(*) 



(20) 



(21) 



Setting E p (£) ^constant (plane wave illumination) , we 
see that in the aberration-free case (</>(xi) = 0) both $ 
and $ raen become delta functions, so that eqs. [T7] and 
[201 both lead to images of the form |G(Afxi)| 2 , where 

M = +1 for eq. [T7]and M = - for eq. HOI 

The advantage of the proposed setup is clear in the 
case where the aberrations are all of odd order: in this 
case, the factor of $ in eq. [20] distorts the image, whereas 
the image formed via eq. [T7]is unaffected. However, the 
price to be paid for this is seen by considering the case 
where the aberrations are entirely even order: the fac- 
tor of 2 in the exponent of eq. [15] doubles the effect of 
the even-order aberrations (compare equation (jT5"j) to eq. 
(|2Tjl ). The reason for this is clear: both the signal and the 
idler contribute to the image and both gain extra phases 
from the aberration. For the odd-order terms, the two 
phases cancel, while in the even-order terms the phases 
add constructively. So, while the setup proposed here 
eliminates odd-order aberrations, it does so at the ex- 
pense of worsening even-order aberrations. This method 
will therefore be of maximal benefit when the aberration 
term of greatest importance is of odd order and is much 
larger than any of the even order terms. 



A further advantage is the improvement in sensitivity 
over standard single-detector imaging due to reduced ef- 
fects of noise. For example, we can look at the effects 
of the detector dark current on the coincidence rate. We 
assume (i) that all fluctuations have zero mean, and (ii) 
the dark current fluctuations in one detector are uncor- 
related with the signal and dark currents in the other 
detector. Then we can write the total detected current 
as the sum of the signal and dark currents (I s and Id), 
and (because of assumption (i)) further split each signal 
and dark current into mean and fluctuation parts: 



Ij = Isj + Idj = (Isj) + (Idj) + SI s j + SI d 



(22) 



where j = 1,2 labels the detector. If the average back- 
ground counting rate is subtracted off, the remaining co- 
incidence rate is proportional to the correlation function 
between the measured currents in the two detectors, 



G^(I u I 2 ) = (hh)-(h)(I 2 ). 



(23) 



Substituting eq. [55]into eq. [23]and making use of the two 
assumptions above, a few lines of algebra quickly shows 
that the effect of the dark current completely cancels out 
of the correlation function. 

The effect of quantum noise is more complicated, and 
has been studied recently in a number of papers [22h25| . 
Since the effect of quantum noise is strongly dependent 
on the parameters of the experiment, it is hard to draw 
general conclusions about its relative effect on quantum 
versus classical imaging. In some cases quantum imaging 
offers benefits on this front as well, but the situation has 
to be examined on a case by case basis. See the references 
cited above for more detail. 



IV. ABERRATION-CANCELLATION WITH A 
CLASSICAL SOURCE 

The odd-order aberration cancellation effect in the pre- 
vious section occurs because the beams strike the lens in 
a spatially anticorrelated manner, leading to the struc- 
ture 0(x') + 4>(x") = <p(x') + </>(— x') in the exponents. 
The entanglement of the beams plays no role here. The 
same effect may be produced by any method that re- 
quires light to strike diametrically opposite points on the 
lens simultaneously. 

To illustrate this, we display in fig. [6] an apparatus 
with classical illumination that achieves the same effect. 
A narrow beam illuminates a rotating mirror, which re- 
flects the beam onto a beam splitter. (The rotating 
mirror could be replaced by any form of beam-steering 
modulator.) The reflected and transmitted beams leave 
the beam splitter with opposite transverse momenta ±q, 
striking the lens on opposite sides of the axis. Over time, 
as the mirror rotates, the entire area of the lens is filled. 
As before, the object is placed in front of the bucket de- 
tector Z?2 ■ The impulse response function for the full sys- 
tem is the product of those of the two individual beams, 
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ration in a simple optical imaging system. Specifically, 
the scheme described here will eliminate effects of odd- 
order aberrations induced by the optical system at the 
expense of amplifying the effects of even-order aberra- 
tions. The results here are in a sense complementary 
to those of [uj - fl5l |. which involved cancellation of even- 
order aberrations induced by the object itself (not by the 
optical system) in interferometry by means of entangled 
photons. It remains for future investigations to see what 
additional types of manipulation of phase effects may be 
possible using variations on the methods considered here. 



FIG. 6: (Color online) Implementation of scheme of fig. 2 
using narrow classically correlated beams for illumination. 



once again leading to the same </>(x') + (/>(— x') structure, 
and to odd-order aberration cancellation. 



V. CONCLUSIONS 

In this paper, we have shown that correlated photon 
methods may be used to manipulate the effects of aber- 
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